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Using the analytic properties of two-point functions in QCD, as well as unitarity, 
bounds on the B meson form factor F(q 2 ) can be derived. Heavy quark spin sym- 
metry, correctly taken into account, is shown to improve these bounds significantly. 



1 Introduction 

An evaluation of the Cabibo-Kobayashi-Maskawa matrix element V c b can be 
obtained by studying the weak decays of B mesons to D mesons: B — ► Dlv. 
The amplitude for this process is: 



„2 _ TUT 2 

q lvl w 



(1) 



xV cb < D(p')\c^ y(l - l5 )b\B(p) > 

The corresponding total cross-section can be taken from experiment. The 
leptonic vertex can be computed perturbatively; to find V c b all there is yet to 
do is to find a way to evaluate the hadronic matrix element. 

We can do so by taking advantage of the fact that B and D are heavy 
mesons. It was shownld that the form-factors which describe such transitions 
are, in the infinite mass limit, expressed through one univer|Sal function: the 
Isgur-Wise function. Using the heavy quark effective theory El (HQET), finite 
mass corrections suited to each particular case can be computed. 

In this paper, we will concern ourselves with the following matrix element: 



< B(p')\brb\B(p) >= (p + p'YF(q 2 ) q = p-p' 



(2) 



Specifically, we will try to obtain bounds on the slope (charge radius) of 
the relevant form factor F near zero momentum transfer : 



2 *o 
P =2 



dF 



,*0 = 4m B 



(3) 



Thus, information on the Isgur-Wise function can be obtained; actually, to a 
good approximation F coincides with this function. 



1 



The method applied in the following makes use of the analyticity and 
unitarity of the theory. A dispersion relation for the vacuum polarization 
function it is written. A Lehmann spectral reprezentation can be used for 
the imaginary part of this function. Keeping in the spectral sum only the 
contribution of certain states we get an integral inequality which relates the 
form-factors along the unitarity cut (physical region for pair production) and 
the value of the function far from this region. The spin symmetry results of 
HQET are used in order to get a better inequality. Finally, applying standard 
techniques related to vector- valued analytic functions!!] the required constraints 
on (||) are derived. 

In the next section we outline the method, using as an example the simplest 
case (taking into account only the contribution of B — B states in the spectral 
sum). In the following one, we address the problem of including the B — B* and 
B* — B* states; for this, the results of HQET have to be correctly taken into 
account. Finally, wc present some numerical results and a short discussion. 

2 Description of the method 

Let's consider the vacuum polarization tensor: 

W{x) = i < 0\T{V(x)V u (0)}\0 >, V" = h"b 
Going into momentum space, we get the invariant amplitude Tl(q 2 ) : 

n^( g ) = (gv-<ftr)n(« a ) 

The analytic properties of this function (in the complex q 2 plane) are easy to 
derive. In the spectral sum 

Imll(q 2 ) = ^ir(2ir) 3 x (4) 

$3<s(q-Pr)WI - E r) < o|v(o)|r >< r|v^(o)|o > 
r 

only states which contain a b and a b quark will contribute. 

Such states can be uniparticle states: Ti, T2, T3, T4 (bound states of the 
b — b pair) or multiparticle states : a B — B,B — B* or B* — B* pair (a B is 
a bound state of the heavy quark b and one of the light quarks u, d or s; the 
difference between B and B* is in their internal quantum numbers; thus B is 
a pseudoscalar particle, while B* is a vectorial one). 

From the above relation, as well as the masses of these particles, it can 
be seen that the function II will have three poles in q 2 = m%< , m^- and ra 2 ^. 
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and a cut starting from q 2 = to = Am? B ( rrij. > to, so this pole is covered by 
the cut). Therefore, we can write a dispersion relation for the derivative of II : 



dq 2 7T J m 2 _ t (t - q 2 ) 2 

1 1 

for g 2 on the real axis and q 2 < rrij. . 

Furthermore, it can be shown that the contribution of a certain type of 
particles in the spectral sum is positive; so , keeping only the contribution of 
B — B states in (^) we can write an inequality like: 

Imll(q 2 ) > ^ 7r(27r) 3 x (6) 
]T 6(q-p T )6(\q°\-E r ) <0\V^\BB><BB\V^\0> 

T(BB) 

The matrix elements appearing in this inequality can be expressed with the 
form factor defined in (||); thus we get : 

ImU(q 2 = t )>^-(l- t -^j ^ \F(t)\ 2 9(t to) (7) 

where nt = 3 comes from the fact that we have three different types of B 
mesons (we presume that the form factor is the same for B u , Bd and B s ) ; 
inserting in (|J) : 



nv)>^/ o ^^mrdt (8) 

In the above relation, there is only one unknown quantity : the F form 
factor (H'(q 2 ) can be computed perturbatively for q 2 <C rrij, ). To get infor- 
mation from (|^) on the charge radius (|^) , we apply the following procedure: 

Defining 

we perform a conformal mapping: 



z= VJ^t-Vu> (10) 

v^o — t + v^o 
3 



which brings the t — q 2 plane inside the unit circle; the threshold to goes into 
z = — 1, the points on the physical cut go onto the circumference of the circle 
; (g) becomes 



with z = e lS and 



1 



^ )= 16^ (1+ ^ (12) 

in the particular case when q 2 = 0. 

Second, we get rid of the singularities of the integrand in the unit disk. It 
can be shown □ that the form factor F also has poles at the square masses of 
the three T particles; in the z plane, they will appear as simple poles z\,z%, Z3 
somewhere on the real axis between z = — 1 and z = 0. To make them 
disappear, multiply the integrand in ( |Tl| ) by the so-called Blaschke functions: 

B(z) = f[^L (13) 
t\ 1 " zz l 

which have zeroes in zi,z 2 ,z 3 , and, moreover, = 1 on the unit circle 

\z\ = 1. Then : 

/■2tt 

1= mz)F{z) B(z)\ 2 (14) 

Jo 2?r 

the integrand being an analytic function. 

Finally, let's consider the following expression : 

1=^- r r \G{z)\ 2 d6, z = e ld . 
27r Jo 

If G is analytic inside the unit disk we can expand it in a power series: 

G(z = e 10 ) = c + Cl e l9 + c 2 e 2 ^ + . . . 
and actually perform the integration : 

/ = M 2 + | Cl | 2 + |c 2 | 2 + ... 



a if we'd knew the residues of F in its poles we would be able to use abetter method - in 
essence, substract the singularities, instead of multiplying them out; seeEl 
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so that 1> I implies : 

l>|c | 2 + | Cl | 2 = |G(0)| 2 + |G'(0)| 2 

This is a Schur-Caratheodory type inequality; in our case it reads: 

1 > |B$F| 2 (0) + |(S$F)f(0) (15) 

Using the normalization F(0) = 1, this quadratic inequality will give a superior 
and an inferior bound on p 2 . 

3 Heavy quark spin symmetry 

The next step is to try to improve these bounds. Obviously, one way of achiev- 
ing this is to take into account as many terms as possible in the right hand 
side of (Q). It is easily seen that, with some phenomenological input,, the 
contribution due to uniparticle states can be computed; thus we getQ 

n fa ) ^ w 22 {q 2 _ + ifa? y to t3/2(t _ g2)2 I^WI ^ ( 16 ) 

where the widths Ty. are physically measurable quantities defined by0: 



cr(e+e~ -> Ti) = 127r 2 (5(t - ) 



Mr 



Further, we try to include contributions from the B — B* , B* — B and 
_B* — B* states. The relevant matrix elements in (||) can be expressed through 
the following form-factors : 

< B*tf t e)\V»\B(p) >= e v p' a ppV(q 2 ) (17) 

tub + m B 



< B*(p', e>)\V»\B*(p, e) >= F x ( g 2 )(e ■ e')P, + F 2 (g 2 )[e M (e' ■ P) + e^e • P)] 

+ F 3 ( q 2 ) ie - P) i e '- P) P, + F,(q 2 )[e a (e> ■ P) - e' a (e ■ P)] 9 " V = t£ (18) 

(where P = p + p', q = p — p'); thus, instead of (P),(^6|) we have : 

27 _ M r IV 
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t 3 / 2 (t-q 2 ) 2 



where ig = (to# + tub*) 2 ,t$* = Am 2 B , ; is a linear combination of F\, 



F 2 and F 3 (see (|24|) below). [| 

Having six unknown functions in it, this inequality is, in this form, of no 
use. At this point, the results of HQET can be of help; it can be shown that 
the form-factors are related in the neighborhood of zero- recoil point as follows: 

V(t),F 2 (t)^F(t),F 1 (t)^-F(t) (20) 

F 3 ,F 4 (t) -> , when t = q 2 -> 

If one assumes i that these relations hold on the entire unitarity cut then 
the inequality ( |l9| ) can be written in terms of a single form factor F. But this 
assumption was shown not to hold, especially near thresholds 

Treating each term in ([l9]) separately, we can still obtain an inequality like 
@: 

1 > |S$F| 2 (0) + \(B<f>F)f(0) + \B$ v V\ 2 (0) + \{B$ v V)f(0) + ... (21) 

The full expressions for <I>, are given in 111 . This relation contains only 

the values of the form factors and of their derivatives at q 2 = ; therefore the 
relations (|2(]) hold, and (|l9|) will become : 

a(p 2 ) 2 - 2bp 2 + c < 1 (22) 
a = 64^(5$) 2 (0) b = 8^(B$)(0)(£$)'(0) 

c = Y,[(.m 2 (o) + my 2 (o)] -1 

with solution : 

b - \/b 2 - ac < p 2 <b + \/b 2 - ac (23) 

This is our final result. 

In deriving ( |l9[ ) , there is a technical point worth discussing. The contribu- 
tion of B — B and B — B* in the spectral sum ([|) are parametrized by single 

b equation rJiih here differs from eq. (13) in by a factor 2 in the term containing V(t), 
because in til the contribution of B* — B states was not included 
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form factors F and V; which appear in the integral @ as \F\ 2 and |V| 2 ; but 
the matrix element of B* — B* is parametrized by four form factors, which 
appear in ( |l9| ) as a positive definite quadratic form. To be able to apply the 
method described in section 2, we have to write this quadratic form as a sum 
of squares; it actually appears this way in (^), F3 being defined as : 

9i 9t It 9t 

W) = (=r - + J^F 2 (t) + — (— - l)F 3 (t) (24) 

If) If) If) if) 



4 Results and discussion 

To obtain numerical results for the bounds (|2^) there are two parameters whose 
value has to be chosen (between certain limits). The first one is q 2 in ([l9|). 
Theoretically, we can chose any value for q 2 from —00 to m 2 ^ ; the best results 

are obtained with q 2 as big as possible. In practice, II' (q 2 ) is evaluated!'! (up 
to three loops) perturbatively; the reliability of this evaluation increases if q 2 
is far from the physical region; it seems that q 2 — 50 GeV 2 is as close as we 
can get to m 2 ^ = 89.5 GeV 2 and still believe in the perturbation series. The 
second parameter is the mass of the b quark, which appears in the evaluation 
of II. As this quantity is not well defined, we choose to vary it from 4.7 to 

5 GeV. Higher mass gives stronger limits. 

In previous work, taking into account only the contribution of T and B-jJ3 
states (not using spin symmetry at all) the following limits were obtainedlJ : 
—5.0 < p 2 < 4.5. Taking into account the contribution oiB — B and B — B* the 
bounds become —0.90 < p 2 < 2.60ll3.These results were obtained with q 2 = 0, 
rn b = 5 GeV and using one loop approximation to compute IT (0). Because the 
two and three loop contribution to II' (0) increase its value by approximatively 
40%E!J (!), these limits are actually stronger than what they should be. 

Using the full apparatus presented in this paper, we obtain : at q 2 = 0: 
-0.2 < p 2 < 1.85 for m b = 4.7 GeV; -0.1 < p 2 < 1.76 for m b = 5.0 GeV; 
at q 2 = 50 GeV 2 : -0.0 < p 2 < 1.6 ,m b = 4.7 GeV and 0.3 < p 2 < 1.2 for 
m b = 5.0 GeV. 

These results show that using only general properties of QCD (as analytic- 
ity and unitarity) we can derive nontrivial constraints on the behavior of heavy 
mesons form factors at transfer momentum close to zero. The use of heavy 
quark spin symmetry brings significant improvements. 

Further results are expected from the application of this method-to_cither 
transitions, like the physically interesting case of B to D* transitionli3,lii 
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